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We develop a non-perturbative "gauge Mean Field Theory" (gMFT) method to study a general 
eflFective spin-1/2 model for magnetism in rare earth pyrochlores. gMFT is based on a novel exact 
slave-particle formulation, and matches both the perturbative regime near the classical spin ice limit 
and the semiclassical approximation far from it. We show that the full phase diagram contains two 
exotic phases: a quantum spin liquid and a coulombic ferromagnet, both of which support deconfined 
spinon excitations and emergent quantum electrodynamics. Phenomenological properties of these 
phases are discussed. 



Amongst the celebrated exotic phases of matter, of 
particular recent interest are the Quantum Spin Liq- 
uids (QSLs) [1]. Behind seemingly innocuous defining 
properties -strong spin correlations, the absence of static 
magnetic moments, and unbroken crystalline symmetry-, 
QSLs display the consequences of extreme quantum en- 
tanglement. These include emergent gauge fields and 
fractional excitations, which take these states beyond the 
usual "mean field" paradigm of phases of matter. Not 
only are these phases challenging to predict and describe, 
they have also proven very hard to find in the laboratory, 
rendering their search and discovery even more tantaliz- 
ing. 

A consensual place to look for QSLs is among frus- 
trated magnets [T] . Frustration allows the spins to avoid 
phases where they are either ordered or frozen, with rela- 
tively small fluctuations and correlations between them. 
Recent experiments have given compelling evidence of 
a QSL state in certain two-dimensional organic mate- 
rials [5], but both microscopic and fully consistent phe- 
nomenological theories are lacking. By contrast, classical 
spin liquids have been conclusively seen and microscop- 
ically understood in the spin ice pyrochlores [3]. This 
raises the possibility, suggested experimentally @] and 
theoretically [5], of QSLs in those rare earth pyrochlores 
in which spins are non-classical, supported by recent re- 
sults on Yb2Ti207 Hj. However, for any material, only 
detailed, quantitative theory predicting the type(s) and 
properties of QSLs that appear and matching experi- 
ments can take the physics to the next level. 

We take up this challenge here for quantum rare earth 
pyrochlores. Our analysis confirms that a "t/(l)" QSL 
phase exists in the phase diagram (Fig. [T]) of a spectrum 
of real materials, and is furthermore supplemented by 
another exotic phase, a Coulombic ferromagnet, which 
contains spinous, but displays non-zero magnetization. 
We also study the confinement transitions out of these 
Coulomb phases, which are analogous to "Higgs" transi- 
tions [6 . Finally, we discuss experimental signatures of 
the C/(l) QSL, and of the [/(I) Coulomb ferromagnet. 

The most general nearest-neighbor symmetry-allowed 
exchange Hamiltonian for spin-1/2 spins (real or effec- 
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FIG. 1. Gauge mean field phase diagram obtained for J±± = 
and > 0. "QSL", "CFM", "FM" and "AFM" de- 
note the U{1) Quantum Spin Liquid, Coulomb Ferromagnet, 
standard ferromagnet, and standard antiferromagnet, respec- 
tively. Phase boundaries with/without white lines indicate 
continuous/discontinuous transitions in gMFT. Note that the 
diagram is symmetric in J^± — ^ — Jz±. 



five) on the pyrochlore lattice is 

H = Y, [j^^Sm ~ J±(S+S- + S-S+) 

+ J±± bjS+S++7*,.S-S7] 

-Kj,±[sr(o,s+ + c*S7) + ^^j]], (1) 

where 7 is a 4 x 4 complex unimodular matrix, and 
C = —7*. The explicit expression of 7 and of the local 
bases whose components are used in Eq. ([T]) are given in 
the Supplementary Material. The first term (we assume 
in this paper J^z > 0), taken alone, gives the highly frus- 
trated classical nearest-neighbor spin ice model, which 
exhibits an extensive ground state degeneracy of "two- 
in-two-out" states. 

In fact, this model has been studied theoretically in 
the special case Jz± = J±± = 0, where it reduces to an 
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"XXZ" model with global XY spin-rotation symmetry 
[7]. There, it was shown that for J± <C Jzz, it is per- 
turbatively equivalent, order by order, to a lattice C/(l) 
gauge theory, with gauge fields that describe the spin con- 
figurations constrained to the spin ice manifold of ground 
states. This gauge theory was furthermore argued to ex- 
hibit a so-called "Coulomb phase" , which corresponds to 
a U{1) QSL phase. Subsequent numerical simulations 
[Hlin] verified this prediction. This Coulombic QSL is not 
only magnetically disordered, but also supports several 
exotic excitations: spinous (called magnetic monopolcs in 
the spin ice literature) , dual "electric monopoles" , and an 
emergent photon. This understanding, however, was lim- 
ited to the perturbative regime J± <^ Jzz and considered 
only the XXZ case. Here we develop a non-perturbative 
method to analyze the full Hamiltonian in Eq. Q. 

Non-perturbative theories of QSLs based on "slave 
particles" have been developed and used extensively in 
SU{2) invariant 5 = 1/2 Heisenberg and Hubbard mod- 
els [10] . Generally these approaches work by embedding 
the Hilbert space on each site in some larger "spinon" 
one, with a microscopic gauge symmetry which acts to 
project back to the physical space. QSL phases are found 
when, in a mean field sense, this microscopic gauge sym- 
metry is incompletely broken in the ground state. Here, 
we follow the spirit but not the letter of these approaches, 
by introducing redundant degrees of freedom not for each 
spin but for each tetrahedron of the pyrochlore lattice. 
This new slave particle representation is, like the afore- 
mentioned standard ones, formally exact, but addition- 
ally naturally describes the Coulombic QSL found before 
in the perturbative analysis, when that limit is taken. 
It also has the added advantage that, unlike in stan- 
dard approaches, the gauge fields appear explicitly in the 
slave particle Hamiltonian, rendering the analogy to lat- 
tice gauge theory more direct and transparent. 

By dint of the theory developed in Refs. IH [71 and [51 
we define our slave particles on the centers of the "up" 
and "down" tetrahedra of the pyrochlore lattice, which 
comprise two FCC sublattices (I/H, with rj^ = ±1) of 



sites, denoted with boldface characters r, of a dual dia- 
mond lattice. The sites of the original pyrochlore lattice 
are bonds of the dual lattice. The perturbative analysis 
of Ref. [7] identified the low energy states of H as the spin 
ice ones, supplemented by spinons corresponding to de- 
fect tetrahedra. As mentioned above, this inspires us to 
enlarge the Hilbert space and define "spinon" slave oper- 
ators, which in turn can be seen as particles in a fluctu- 
ating vacuum (the two-in-two-out manifold dear to the 
spin ice community). We consider T-Lbig = T-Lspin ^ Ho, 
where Hspin = (S*jv^i/2 ^^'^ Hilbert space of Eq. Q 
and Hq is the Hilbert space of a field Qr G Z. Qr is de- 
fined on all the sites of the dual diamond lattice and, at 
this stage, is free and unphysical. We further define the 
real and compact operator ip^ to be the canonically con- 
jugate variable to Qr, [Vr,Qr] = i- In Hq, the bosonic 
operators $J = e*'^"' and $r = e"**^' thus act as raising 
and lowering operators, respectively, for Q^. Note that, 
by construction, |$r| — 1- We now take the restriction 
of T-Lbig to the subspace %, in which 



Qr = ??i 



(2) 



where the e^'s are the four nearest-neighbor vectors of 
an = 1 (I) diamond sublattice site. This constraint 
can be viewed as analogous to Gauss' law, where now Qr 
counts the number of spinons. The restriction of Q^., 
and ^\ to H exactly reproduces all matrix elements of 
the original Hspim with the replacements 



(3) 



Here r e I, and s^^i.sl,., act within the 'Hspin subspace 
of T-Lbig- Note especially that, by itself, s^^., ^ 5^^., is not 
the physical spin, and does not remain within %. 

In this paper we focus on the case where J±± — 
(which otherwise introduces additional complications to 
be dealt with in a separate publication), and the Hamil- 
tonian then becomes 



r6l,II [ rel p,,v^p. 



+ 



V V $ s+ s- 

/ y / y ^r-e,, r-e,^^r,r-e^^r,r-c 
r^II fi.iy^fi 

r^II fijiy^fj. 



(4) 



const. 



The integer-valued constraint in Eq. ^ commutes with 
H and thereby ensures that Eq. (4]) is a U{1) gauge the- 
ory. Explicitly, it is invariant under the transformations 




-Xr) 



(5) 



with arbitrary Xr- This invariance, and the Gauss' law in 
Eq. ^ can be made formally identical to that in lattice 
electrodynamics by writing s^^., = E^r' and s^^., = e^^^'--' , 
where E and A are lattice electric and magnetic fields [7] . 
This clarifies that s^^., is to be regarded as an element of 
the U{1) gauge group. However, the notation is unnec- 
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essary and we use it only when conceptually valuable. 

Eq. g can be viewed as spinous hopping in the back- 
ground of fluctuating gauge fields, and thereby lends 
itself to the application of standard mean field theory 
methods for lattice gauge models [11], which we call 
gauge Mean Field Theory (gMFT). Upon performing 
gMFT, we will get a Hamiltonian for spinons hopping 
in a fixed background. Specifically, we perform the re- 
placement: 



$1'$(s)(s) 



(6) 



($^$)(s)s-2($1'$)(s) (s) 



and thereby split the Hamiltonian in a spinon part , 
and a gauge part . Note that unlike conventional 
Curie- Weiss mean field theory, which entirely neglects 
any quantum entanglement, gMFT, while suppressing 
some fluctuations, still allows high correlations and en- 
tanglement. 

The gMFT order parameters are closely analogous to 
those in C/(l) Higgs theory [HJIT^]. A non-zero expecta- 
tion value (s*) 7^ implies the phase of is relatively 
well-defined, i.e. there are small fluctuations of the vector 
potential A. The converse case, (s^) = would indicate 
confinement, but does not occur here. A non-zero scalar 
expectation value, ($) ^ 0, analogous to a Higgs phase, 
indicates spinon condensation and generation of a mass 
for the gauge field, and a conventional, non-exotic state. 
Combined with (s*) 7^ 0, it also implies "XY" magnetic 
order. Conversely, (<i>) = indicates the spinons have 
a gap, and is characteristic of the Coulomb phase. The 
remaining gMFT order parameter, s^, is gauge invariant, 
and thus indicates only the presence ((s^) ^ 0) or absence 
((s^) = 0) of "Ising" magnetic order, i.e. time-reversal 
symmetry breaking. Combining this together, the phases 
in gMFT are summarized in Table |lj We emphasize that 



TABLE I. Order parameters and phases in gMFT. 
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phase 










QSL 





/o 


/o 


CFM 






/o 


FM 









AFM 



despite the fact that ($) does not appear explicitly in the 
decoupling in Eq. (|6]), the gMFT does generally allow for 
Higgs phases where <i> is indeed condensed. As we will 
show below, the Higgs phase appears in a manner similar 
to Bose-Einstein condensation in an ideal Bose gas. 

We now use the following Ansatz, valid when J± > 
(which we assume hereafter). 



■ sm a 



■ cos 9, 



(7) 



where /i = 0, ..,3 and e = (1, 1, —1, —1), which assumes 
translational invariance and fully polarized "spins" s, in 
accord with Eq. ([6]), and is compatible with FM polariza- 
tion along the (global) x axis ((s+) = (s~)). Note that 
Eq. ([7]) shows that the gMFT allows fluctuations of E 
and A, so long as 9 ^ n/2 and 6* 7^ 0, respectively. Defln- 
ing the dot product, and through it the vector notation. 



u V 



i(uH 



we find 



H. 



MF 



2^2^heff',.(r) 



(8) 



4£^ Jz±A cos6' and h^jj = AJz±Ii sin^ 



where h 

2J±l2COs9, and we have defined /i = ^tii^l^r+e,,) (^^o 

summation implied) and I2 = J2i,^p,{^t^T+e^-e^) (a* is 
fixed). These quantities turn out to be independent of 
the diamond bond /x. To treat the spinons, we relax 
the |$r| = 1 constraint to a global one by introducing a 
Lagrange multiplier A via the term A^^. (|$rP — l) in 
a path integral formulation, with free integration over <i> 
and $*. The spinon Lagrangian is 



/•MF 



E 



■G7 



(9) 



where A^u.c. is the number of unit cells, [Gk^ujn]ab ~ 
(<i>^<i>a), and we find the equal time Green's function to 
be 



G 



k,T=0 



Mk 
|Mk|' 

zt 



(10) 



± 



where Mu = Eu ^) = 

^ ^k ^k 

^k(^.^) = v^A-^±(0), ll{9) = ij±cos2 0LkT 
|iJz±sin26'Mk|, Lk = I],,,i/<^ cos [k • (e^ - Gi,)]. A 
couple of remarks are in order: (i) A > maxk 
(ii) the spinon dispersion relations are a;^(0,A) = 

^m7.z^{9,x). 

The gMFT consistency conditions on 9 and A (for 
fixed J±,Jz±) arise from requiring Eqs. ^ and (sp = 

h^/(2|h^|), and from the normalization condition on 
and can be written 




2J,±/i(0,A) 



2J,±h{9,X)tan9 + J±l2i9,X) 
A) = l 



(11) 



where I3 {9, A) = The explicit expressions for the 

Ii = iV~(^ Ek-^k' needed to solve Eqs. (Ill, are readily 
derived from Eq. ( |10[), and are given in the Supplemen- 
tary Material, Eq. (|B5|). Since Eqs. (11) may allow sev- 



eral distinct solutions, we must choose the solution of 



Eq. (Ill with the lowest energy. In the mean field ap- 



proximation, the ground state energy can be calculated 
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by taking the expectation value of the Hamiltonian. We 
find, per unit cell, ecs = ^av + efem, with 

eav = -2/2(6', A) cos^ 0J± - 4/1 (61, A) sin 26* J,± (12) 
ek^n = l|JoJ+ie,\)+uJ^{9,X)], (13) 

where — y/2J^z^. Here Ckin measures the "kinetic" 
energy associated with the spinon modes, while €av repre- 
sents the "background" energy in which the latter evolve. 

We now discuss how the different phases are obtained 
from the solutions of the gMFT equations. Condensed 
and uncondensed phases are distinguished by the value 
of A. As in the theory of superfluidity, condensation 
is synonymous with off-diagonal long-range order, i.e. 
(<I>r)*($r/) = lim|r_r'|-).oo('i'J<I'r') ^ 0. This expectation 
value ($r'*J) = ^ul. J^k'^kS^^'^"'^"^ is non-zero in the 
long-distance limit if and only if the usual conversion of 
the sum to an integral fails, i.e. if there exists one ko 
such that Gko — 0{Nu.c.)- Like the chemical potential 
in Bose-Einstein condensation, A in a Higgs phase differs 
from its minimum allowed value by a sub-extensive part. 
A = Ainin(^) + ]^— therefore defines condensation, where 

S is of order 0(1), and Aniin(6') — maxi^£^{6). If instead 
A — Amin(6') is 0(1), one has a phase with deconfined 
gapped spinous and a Coulombic gauge structure. As we 
already mentioned above, these classes of phases can be 
further subdivided into "polarized" (i.e. with magneti- 
zation along the local z axis) and "unpolarized" when 
6' 7^ and = 0, respectively. 

The phase diagram resulting from the gMFT solution 
(see Supplementary material) is shown in Fig. [I] It 
contains two "exotic" phases in which spinous are de- 
confined and uncondensed, indicated as QSL and CFM. 
The QSL state, with 6* = 0, is completely absent mag- 
netic order, and is the phase studied in Refs. T' and '8', 
Its low energy physics mimics quantum electrodynamics, 
and thereby contains a photonic excitation (gapless and 
linear near the origin) and gapped fractional monopole 
excitations (spinon and "electric" monopole) that inter- 
act via Coulomb interactions. In the present formal- 
ism, the photon is only obtained once quadratic fluctu- 
ations around the gMFT solution are considered, but is 
a universal feature of the exotic phases. The CFM, or 



"Coulombic Ferromagnet" phase, is a new phase of mat- 
ter that can be seen as a polarized version of the C/(l) 
QSL. Despite being magnetic, its elementary magnetic 
excitations are spinous rather than spin waves, and it 
also supports a gapless photon mode. Indeed, in gMFT 
the transition from the QSL to CFM is second order, and 
consists of a continuous rise of magnetization from zero. 
For larger Jz±,J±, one obtains Higgs phases, which are 
conventional states of matter without exotic excitations 
and are continuously connected to the usual magnetically 
ordered states described by Curie- Weiss MFT. Interest- 
ingly, we find the exotic CFM state is considerably more 
stable than the "pure" QSL, occupying a much more sub- 
stantial portion of the phase diagram. 

How do we recognize a Coulomb phase in experiment? 
A generic sign of fractionalization is a two-particle con- 
tinuum in inelastic neutron scattering, two spinous being 
excited by one neutron [U 113) . In addition the photon 
can be detected directly by inelastic neutron scattering, 
as a linearly dispersing transverse mode. It is, how- 
ever, more challenging to observe than the usual acous- 
tic spin wave, because its scattering intensity becomes 
small (oc w) at low energy (see Supp. Mat.), in con- 
trast to the spin wave for which the intensity diverges 
(^ in the same limit. Interestingly, the pinch points 
in the static structure factor present for classical spin ice 
are absent for the quantum Coulomb phase [S], so this 
is not a useful measurement. Perhaps the most strik- 
ing signature of the Coulomb phase is likely to be ther- 
modynamic. Like the phonons, the photons contribute 
as BT^ to the specific heat at low temperatures, but 
their speed is Vphoton ^ J <^ c, the speed of sound. 
Crudely estimating J ~ 2 K appropriate for Yb2Ti207, 
we obtain a coefficient Bphoton ~ lO'^mJ/mole-K^, ap- 
proximately 1000 times larger (!) than the phonon con- 
tribution Bphonon ~ 0.5 mJ/mole-K^ measured for the 
isostructural material Y2Ti207 [14) . 

With a phase diagram and a new phase of matter in 
hand, we take heart at discovering yet more new exciting 
facts in the pyrochore lattice. Future studies should ad- 
dress the more frustrated case J± < 0, phase transitions 
in applied field, and the influence of defects. 

We thank Peter Holdsworth, Bruce Gaulin and Kate 
Ross for discussions. This work was supported by the 
DOE through BES grant DE-FG02-08ER46524. 
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SUPPLEMENTARY MATERIAL 

Appendix A: Geometry 

1. Local bases 

The local cubic bases in which the Hamiltonian Eq. ([T]) 
is expressed are the following (a^ , , ) bases 



eo = (l,l,l)/\/3 
ei = (l,-l,-l)/\/3 
62 = (-1,1,-1)/V3 
. 63 = (-1,-1,1)/V3, 



ao = (-2, 1, 1)/V6 
ai = (-2,-1,-1)776 
SL2 = (2,1,-1)/V6 
^ a3 = (2,-l,l)/V6 



(Al) 



hi — ei X SLi, such that spin Si on sublattice i is Si 
S+(a, - ih,)/2 + S7(a, + ih,)/2 + Sfe,. 
The 4x4 matrix 7 introduced in Eq. ^ is 



7 



f 1 w 

1 w"^ w 

w w"^ 1 
.,2 



(A2) 



where w 



w 1 / 
g27ri/3 jg ^ third root of unity. 



the first Brillouin zone can be considered the "cube" with 
unit sides described by —1/2 < qi < 1/2 (note that the 
Qi's are dimensionless) . 



Appendix B: Brillouin Zone Sums 



In the main text, we defined 

/i = e^($r<f>r+e„) (no summation implied), (Bl) 



l2 = Y,('^l'^r+e,-eJ (//is fixed), 
For brevity, we define the I^'s through 



nr., -t^ii.r. 



(B2) 
(B3) 

(B4) 



keBZ 



q^Ecubc 



(q is defined in Eq. (A3|), whose explicit expressions are 

Iq(^,A) = 



1 fj^ 



2 V 2 [Afqi 



(— 1 — cos 2TTqi + cos 27rg2 + cos iirq^) , 



-A) = 



(cos 2nqi + cos 27rg2 + cos 27rg3) . 



(B5) 



where we used the Fourier transform convention $r 
jv^ Ek *k e'^ '', and in turn 



(<f|:$r'> - ^ 5][Gk]2ie^'"(""'-'^ for r e I,r' e II, 

(B6) 



with Gk defined in Eq. ( 10 1 



These Brillouin zone sums need to be evaluated. While 
in the Coulomb phases, they can be directly identified 
with the corresponding integrals, the situation is more 
complex within the condensed phases. We address this 
issue below. 



2. Lattice vectors 

The four nearest-neighbor vectors of a I-sublattice di- 
amond site (sublattice I corresponds to "up" tetrahedra) 

are — s^e^, where a is the usual FCC lattice spac- 
ing. The four pyrochlore sites of the "up" tetrahedron 
centered at the origin are located at e^/2, ^J, — 0, .., 3. 

The FCC primitive lattice vectors are A^ = eo — e^, 
i = 1, ..,3, while the reciprocal lattice basis vectors are 
defined as usual by Bi = 27r ^^^^-' and its cyclic permu- 
tations, where v^.c. = Ai • (A2 x A3) is the volume of the 
(real space) unit cell. If the g^'s are defined as 



(A3) 



Appendix C: Brillouin zone sums in the condensed 
phases 

The condensed region occurs if /3(Aniin) < 1- In this 
case, we must allow for the subextensive part of A. Let 



A = A„ 



iV2 



(CI) 



where Nu.c. is the number of unit cells and S > 0. Then 



the sums Eq. ( B4 ) cannot merely be replaced by the cor- 



responding Brillouin zone integrals, but must be written 
as follows: 



1 



/(kmin) 

' II.. r.. 



1 



6 



and 



Nu.c. 



/(k) = 



k#k„ 
1 



E 



/(k) 



(C3) 



V; 



BZ 



E 

k5^k„i, 



where Vbz is the volume of the Brillouin zone. In general, 
we denote this decomposion of an integral as / = /min+-^', 
where /' is the continuous integral part evaluated at A = 
Amin, and /niin is the first term in the right hand side of 



Eq. (C2). In /min, we need only retain the part which is 
non-vanishing as N^ ,,. ~^ oo. 

The value Amin is determined by the condition that 
minq — 0, which corresponds to Amin — maxq£q (0) — 
maxq(i cos^ dJ±Lq + \\Jz± sin26'Afq|). Here, this max- 
imum maxq ^q always occurs for the (dimensionless) 
wavevector of the form Qmin = (9ii92i<Z3) — (0,9,<z) 
(qmin corresponds to ko in the main text). Defining 



y ~ cos 27r(7, 



(C4) 



maXq £q is realized for the wavevector (0, q, q) corre- 



sponding to 

1 - 
-1 



1 ( J2lJz±,t 

8 Ji[J±M 



for J2[J,±,9]<4Ji[J±,9] 

for J2[J.±,e] >4Ji[J±,0] ' 
(C5) 



where we defined 

Ji[J±,0] = I J± CDs'" 9, 
J2[Jz±,0] = ij,±sin20. 



(C6) 
(C7) 



i.e. £^{0) = Ji[J±,9]LuT \J2[Jz±,0]M]i\, and assumed 
J2[Jz±,d] > 0. Then we find 

Amin = (C8) 

f 6 Ji [ J± , 6] + for J2 [Jz± , 0] < 4 Ji [ J± , 0] 

1 4 J2 [Jz± ,0]-2Ji [ J± , 6] for J2 [Jz± , 9] > 4 Ji [ J± , 9] ' 

There is one subtlety. When 9 — 0, — , so that 
there will be a degeneracy of the minimum energy state. 
This changes the condensed contributions. The two cases 
can actually be treated together, provided we formulate 
everything in terms of a new variable p, the "condensed 
density" , rather than S, where 



P{9 > 0) 
Pi9 = 0) 




(C9) 



i.e. the condensed density is larger in the case 6* = by 
a factor of 2. Then, from Eqs. (11) and (B5), we find 

P 



h 



h = p{l + 2y, 



I'2- 



(CIO) 
(Cll) 

(C12) 



When = 0, the minimum is reached at qmin = 0, i.e. 
2/min — 1- Eq. ( |C11 1 then reduces to I[{9 = Q) = since 
Ii {9 — 0) vanishes identically (which reflects the fact that 
the two FCC sublattices are decoupled), and Eq. (C12) 
reduces to I2 — 3p - 



I'2- 



Appendix D: Ground State Energy Calculation 

In some cases, we find multiple solutions of the mean- 
field equations. The physical solution is the one with the 
lowest ground state energy. To find it, we must evaluate 
the energy in the mean field approximation. There are 
different ways to address this. Here we show how to cal- 
culate the ground state energy "directly" , i.e. by taking 
the expectation value of the Hamiltonian. We arrive at 
Eqs. (12 1 and (13) given in the main text. 

We write the Hamiltonian in Eq. ^ as 



H = 



E 

rGl.II 



J. 



(Dl) 



We can treat the gMFT as a variational calculation, 
with the trial wavefunction being the ground state of the 
gMFT Hamiltonian. The expectation value of the second 
term is straightforward, but that of the first is not. We 
need to write the Hamiltonian generalization of the ap- 
proximation of "softening" the j^j = 1 constraint in the 
path integral. 

To do this, we write <I>r = Xr + iyr, and introduce 
canonical momenta Pa;^r:J'y,rj such that [a;r,P2:,r] = * etc. 
Then if we write Hr = Px,r+iPy,n the equivalent approx- 
imation in Hamiltonian form is 



H 



smce 



E 

rel,ll 



Jz 



Ulll, + Xi<^>l^,-l)'>+n[^,s], (D2) 



H - 



1 d 



|$| dm V ^1* 



l$l 



d 



because |$| = 1, as enforced by A(<i>J$r ~ 1) with A > 0. 
One can check that, writing the path integral for the 
Hamiltonian Eq. (D2) and integrating out Hr, one ob- 
tains the action used previously. Note that we can see 
that there should be four harmonic oscillators per unit 
cell, arising from Px,r,Py,r for the two sites r in the di- 
amond basis. This implies in turn that each of the two 
spinon branches uj^ is doubly degenerate. 

Returning to the energy calculation, we simply can 
take the expectation value of Eq. (D2). Note that, be- 
cause we choose ($J$r) = 1, the expectation value of 
the second term in the brackets (with the A coefficient) 
is zero. Hence the energy is 

^GS= '^{^l^r) + {n'^,s\)=Ek^n+Eay. (D4) 

rel.II 
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The second term is 

Eay ^ {n[<P,s]) (D5) 
= Nu.c. [-2/2 cos^ 0J± - 4/1 sm29J,±] , 

i.e. half of the constant term obtained from the replace- 
ment Eq. Q (up to a minus sign) . The "kinetic energy" 
Ekin requires more thought. We can calculate it by the 
path integral. One obtains 



(ntn^) = 



2S{t = 0) - 



J zz 



(D6) 



Each of the two terms in the square brackets is formally 
divergent, but together they give a finite answer. This 
can be seen by writing it all in Fourier space: 



J. 



rel.II 



(D7) 



2Nu.c. f duj- 



Now we can rewrite the Green's function as 



{ni + zi){ni + ziy ^^^^ 



where f2„ — aj„/\/2 J^^, and z± — z^- Changing vari- 
ables from LOn to ri„, we have 

'drj„ nl{z\ + zl)/2 + z\z'^_ 



Ekin — 2A^M.c. \f2j'z 



27r {^^^z\)(Sll^zl) 



(D9) 



We note that this formulation makes it clear that this 
form of the energy is variational. Specifically, even if 
9 is not chosen to satisfy the mean field condition, this 
should give an upper bound to the ground state energy. 
This means that this form of the energy should be mini- 
mized when 9 equals its true value, which is the minimum 
energy mean field solution. So this form could be used, 
if desired, to search for the ground state by just mini- 
mizing this energy as a function of 9. Note that we do 
need to choose A as a function of 9 to enforce the |$| = 1 
condition for this to be valid, though. 



Appendix E: Phase Boundaries 

Most generally, the phase transitions of Figure [l] are 
found numerically as the loci of points where the lowest- 
energy solutions to the gMFT equations have different 
{9, A) characteristics on either side of the transition line. 
It turns out that some of these lines coincide with curves 
which can be found more directly. Here we outline which 
transitions have a direct physical meaning, and how we 
find them. 



We define, and shall use extensively, the dimensionless 
exchange constants 



J^ = 



J ZZ 



(El) 



Finally, the problem being symmetric in Jz± — > — Jz±j 
we will always assume Jz± > 0. In the latter case, all the 
values of 9 of interest are contained in the [0,7r/2[ inter- 
val, so that J2, defined in Eq. (C7), is always positive. 
Hence, | J2I = ^2- 



1. QSL - AFM 

In both the QSL and AFM, 6* = 0. This considerably 
simplifies the problem. Using Eqs. (Ill and (B5), we 



find that the maximum allowed J± within the QSL is 
such that A = 3J±, from which 



J± 



QSL 



■^^z ) max 2 




3 — \E\^ 



(E2) 



follows. The actual QSL-AFM transition happens when 

\ ~~ / max 



2. QSL - CFM 

We find that the QSL-CFM transition is continuous in 
9. We can thus find the transition line as follows. We 



expand the consistency equations Eqs. (11 1 to first order 
in 9, around 6* = 0, 

h{9,\) = 9\Jz±\h{\) + 0{9^) (E3) 
/2(e,A) = /2(O,A) + O(02) (E4) 
h{9,\) = h{Q,\) + 0{9^) (E5) 



where Ii 




(note that Ii has the 



dimension of 1/J), so that Eqs. (11) become 



^±/2(0,A) = 2j2±/i(A) 
/3(0,A) = 1 



(E6) 



where neither /2(0, A), 13(0, A) nor /i(A) depend on Jz±- 
For every value of J±, A is the solution of 13(0, A) = 1, 
so that we find Jz±/ Jzz at the transition to be 



Jzz / crit 



QSL-CFM 



^ J±/2(0,A) 
2J../l(A) ' 



(E7) 



with A(J±) first determined by solving 13(0, A) = 1 at 
fixed J±. 
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CFM - FM 



CFM - AFM and AFM - FM 



The CFM-FM transition coincides with the first ap- 
pearance of a condensed solution (at J± fixed for in- 
creasing Jz±) to the gMFT equations Eqs. (11). To 
find this transition Hne, we use the generic form of 
£^ = JiLk T |>/2-Wk|- In other words, we treat Ji and 



J2 defined in Eqs. (C6) and (C7) as parameters. We also 



define r = J2/J1 = tan^, so that 



(E8) 



The system being at the onset of a possible condensate, 
we know the possible values of A at this onset (see above 
section, Eq. (|C8|), 



'A = 2Ji (2r- 1) 
or A = Ji f 6 + ^ 



so that, from 13(6', A) 
either 



(E9) 



1 (with these values of A), we get 



Jiir) 



, L 



.til v/K4 + i'|Mk|)-2-ikJ ' 

(ElO) 



or 



Jiir) 



u=±i y'6-Lk + r(ir + z^|Mk|) 

(Ell) 

i.e. Ji as a function of r. Using the above equations, we 
can rewrite the consistency equation Eq. (11) in a form 



suitable for determining r as an unknown. For example, 
in the first case Eq. ( |E10[ ), the transition line is best 
found as a function of J±, so we fix J±, and use the 
simple relations (obtained from Eqs. (C6) and (C7|) 



arccos 



2 tan( 




J4 



(E12) 
(E13) 



Using these relations in Eq. ( 11 ) and Eq. (ElO ), we obtain 
an equation which can be solved numerically for r as a 
function of J±. Then, by substituting back, we get Jz± 
and hence the phase boundary in the J±-Jz± plane. In 
the second case Eq. (Ell), we fix 9 (as opposed to J±), 
use the relations 



J4 



2Ji 



cos^ 9 

Jz± = \r--^J±, 
2 tan 6/ 



(E14) 
(E15) 



and proceed similarly. 



The transitions are here discontinuous in 9, and we find 
the transitions indirectly, except for a small portion of 
the CFM-AFM transition line, which coincides with the 

j2 

limiting boundary described above, with A = 6J1 -|- i^j-. 



Appendix F: Benchmarking gMFT 
1. Classical approach 

In Figure [2j we compare the phase diagram of the 
gMFT with that expected by treating the spins in Eq. ([T]) 
classically. The large J±/ Jz± region is found to be anti- 
ferromagnetic, while the large Jz±l J± is ferromagnetic. 

The transition line is found analytically. To achieve 
this, we use the fact that, both in the FM and AFM 
phases, the classical spins take a specific and simple form. 

In the FM phase, 



5™ = (a;ao+ybo-l-zeo)/2 
Sf^ = {-xai-yhi- zei)/2 
SfM = (-xa2-z/b2-ze2)/2 
(a;a3 + yh^ + ze3)/2 



(Fl) 



cFM 
•^3 



where x, y, z are free parameters with + + = 
1. In polar coordinates, we write x = sin0cos0, y — 
sinGsinc/), z = cos0, and the energy Eq. ([T]) per unit 
cell becomes 



FM 
^class 



1 



-2cos^e.hz+4:Sm^ej± (F2) 
+4 sin 20 ^cos + sin Jz± . 



The energy is then found to be minimized for (p — tt/3 
and O — arctan 



l+2j±-.^H-4j±(J±-|-l)+64j2 



which 



yields e™ , 



Jz 



-1 + 2J± - a/(1 + 2J± 



where Ji = Ji/Jzz, Eq. (|Elj). 

In the AFM phase, the solution takes the form 



Yh, 



(F3) 



with arbitrary X, Y such that X'^ +Y'^ — 1. This gives 
the energy per unit cell 



,AFM 



3J± 



(F4) 



Equating e™j,^ = e^™, we readily find that the transi- 
tion line takes the form 



z± 



class 



'J±(6J±-1) 



(F5) 
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The comparison between the classical and the gMFT di- 
agrams shows that the two approaches agree in the semi- 
classical limit of the gMFT, i.e. for J±l Jzz-, Jz±/ Jzz ^ 1- 



2. Pert ur bat ive approach 

In the small Jz±/ Jzz find J±/ Jzz limit, the perturba- 
tive approach developed in Refs. [7] and S] applies. There, 
the effective Hamiltonian is 



rrport 



{i j,/c,/,m,n}— O 



h.c. 



-Ji 



(3) 



with K 



12Ji 

— ^ and J(3) 



J 



(F6) 



and O represents the 



(flat) hexagons of the pyrochlore lattice. While the first 
term induces the spin liquid physics, the second, taken 
alone has the six ferromagnetic ground states described 
in Ref. i4, with polarization along the (100) directions. 
The actual phase diagram is therefore expected to differ 
from that of the gMFT in the small Jz±/ Jzz region (but 
agree when Jz±/ Jzz = 0). Dimensionally we find that the 
transition between the FM and the CFM in this region 
occurs along 



z± 



PT 



zz / crit 



J± 

J ZZ 



3/2 



(F7) 



in the region where perturbation theory applies. A sketch 
of the expected gMFT diagram which includes this per- 
turbative limit is sketched in Figure |3] 



Appendix G: Physical Properties 
1. Staggered magnetization 



where r and r' are on diamond sublattice I and, of course, 
i e (r, r + e^) and j € (r', r' + e^), and where we made 
use of Wick's theorem. From Eq. (10), we find 



1 



^Z+(0,A)e*'' ("-'-'-). (G3) 



In the QSL, the sum can be replaced by its corresponding 
integral, 

-^5]z+(0,A)e''' ('-'-"-)-> / Z+(0,A)e''' ('-'-"-)=0, 

(G4) 

because |r' — r| — > oo, and never reaches within 
the integration domain. We thereby recover (S^) = 
in the QSL. In the AFM, however, the "minimum term" 



survives via the subextensive part of A, Eq. (|C1|, so that 

($!$w) : 



J zz^ 



(G5) 



since k,„in = when 9^0. Using Ec^. (CIO), we find 
^1 



2 S 



1 



= 1 




as defined in Eq. (|E2]). Similarly, ((i>r+e„*J/_ 



^ I in the AFM, so that 




(G8) 



in the AFM phase for Jz± — 0, where we have defined 
J± = J±/Jzz and J^ = {J±lJzz)^^^- (S^) at J,± = is 
plotted versus J± / Jzz in Figure H] 



Here we calculate the "staggered" magnetization 
within the QSL and AFM phases, i.e. the local mag- 
netization of each spin. This quantity will allow us to 
distinguish between the QSL (uncondensed) and AFM 
(condensed) since the reasons for overall zero magnetiza- 
tion in both phases are different. Indeed, in the AFM, 
the local the spin expectation values are non-zero but 
compensate, while in the QSL, the (local) XY symmetry 
remains unbroken, and (5^) = everywhere. 

We calculate the staggered magnetization along the lo- 
cal X axes. To do so, we compute the correlation function 
(S,^S~), where sites i and j are taken far apart from one 
another, where they can be considered independent. 

(S + ST) = (<f]:s^r+e,*r+e,$|:,+e,.Sr;,r'+e„<fr') (Gl) 
= (s+)(s;)(<l)t$rO(<i>r+e,<i>^+eJ, (G2) 



2. Structure factor 

The excitations above the uncondensed ground states 
of Eq. (|4| comprise the photon, spinous and "electric" 
monopoles jHIT]. In general, the spinous and monopoles 
are gapped excitations, and so do not contribute to the 
low-energy part of the structure factor. Below some 
threshold, only the photon survives, and gives rise to 
a universal form of the scattering. Here we sketch the 
general argument which shows that it appears in the in- 
elastic structure factor as a linearly dispersing mode at 
arbitrarily low energy, with an unusually small spectral 
weight and polarization dependence which distinguishes 
it from the usual spin waves. We will be deliberately 
vague here about lattice details, orientation of local spin 
axes, etc, as these affect only 0(1) dependencies of the 
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FIG. 3. Modified gMFT diagram, which takes into account 
the known perturbative limit J^^jJ^^ <^ 1. Note that the 
FM-CFM transition (white dashed hne) in the latter region 
is a sketch. 



amplitudes and not the important scaling properties. A 
much more complete and microscopic study of the struc- 
ture factor in the QSL and CFM phases is planned for a 
future publication. 

Rather than derive it microscopically, we rely on the 
universality of the low energy description of the emer- 
gent photon. We expect that it is described by the usual 
electromagnetic action, which may be written in the con- 
tinuum for small momentum and low energy. We write 
this in Euclidean time in terms of the vector potential, 
and choose the gauge in which the scalar potential van- 
ishes. The result is familiar: 



5, 



QED 



Cl 



{drA) 



C2 



(V X A) 



(G9) 



Here ci and C2 are phenomenological parameters related 
to the effective dielectric constant and magnetic perme- 
ability for the effective electrodynamics. Importantly, the 
photon velocity is given hy v = \fc2fc\. We recall the 
relation between the formalism of Eq. Q and that of 
lattice electrodynamics, given in the main text. 




0.2 0.4 0.6 0, 



1.0 1.2 1.4 1 



FIG. 4. Staggered magnetization (S^) as a function of J±/ Jz 
Eq. ([G8|, for J^± = 0. 



Srr' - Kt' - Err' 



(GIO) 



where r is on diamond sublattice I. Because the trans- 
verse components of the spins, S*, involve the spinon 
fields, we expect that they contribute only above the two- 
spinon threshold. For low energies, it suffices to consider 
only. We take the continuum limit by defining the 
electric field to be "oriented" along the diamond lattice 
bond I — II on which the spin lies and centered midway 
along it, i.e. 



E. 



+e„/2 



E, 



■r,r+ej,c^. 



(Gil) 



Similarly, recall that the magnetic moment due to 
is along the direction in real space, 



(G12) 



So we see that, in this low energy subspace where excited 
spinous may be neglected, the spin operator is precisely 
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the electric field operator. Then, using the usual relation 
E(r) = -VrV^(r) - drA{r) for y = (gauge choice), 



(G13) 



where aj„ is the bosonic Matsubara frequency. Thus the 
Matsubara correlation function of the spins is 

{SU.^o.Ji.J^^l{^-K-^Al.J, (G14) 

Calculation of the gauge field propagator is a textbook 
exercise: 

(^-.-.AL) = cr^%^^. (G15) 



Inserting this into Eq. (G14) and analytically continuing 
iun — > w + iO~^ , we arrive finally at the inelastic structure 
factor 



J^is.,uj = - Im {S- 



o+iO+ 



(G16) 



ci 



uj S{uj — t;|k|). 



We see that the photon appears as a sharp peak in the 
structure factor, with a weight proportional to its fre- 
quency u). This indicates, as mentioned in the text, a 
strong suppression of the weight at low energy, especially 
relative to the familiar case of a spin wave, for which the 
weight diverges like 1/uj in the same limit. 



